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ABSTRACT 

Spectral and photometric observations of nearby galaxies show a correlation 
between the strength of their mid-IR aromatic features and their metal abun- 
dance, and a deficiency of these features in low-metallicity galaxies. The aro- 
matic features are most commonly attributed to emission from PAH molecules. 

In this paper, we suggest that the observed correlation represents a trend of PAH 
abundance with galactic age, reflecting the delayed injection of PAHs and car- 
bon dust into the ISM, by AGB stars in their final, post-AGB phase of their 
evolution. These AGB stars are the primary sources of PAHs and carbon dust 
in galaxies, and recycle their ejecta back to the interstellar medium only after a 
few hundred million years of evolution on the main sequence. In contrast, more 
massive stars that explode as Type II supernovae inject their metals and dust 
almost instantaneously after their formation. After determining the PAH abun- 
dances in 35 nearby galaxies, we use a chemical evolution model to show that 
the delayed injection of carbon dust by AGB stars provides a natural explana- 
tion to the dependence of the PAH content in galaxies with metallicity. We also 
show that larger dust particles giving rise to the far-IR emission follow a distinct 
evolutionary trend closely related to the injection of dust by massive stars into 
the ISM. 
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1. INTRODUCTION 

Spectral and photometric observations of nearby galaxies with the Infrared Space Obser- 
vatory (ISO) and the Spitzer Space Telescope have provided the opportunity to investigate 
the inter-relations between global galactic properties, such as morphology, star formation 
activity, spectral energy distribution (SED), metallicity, and dust abundance and compo- 
sition. In particular, these observations have enabled detailed studies of the correlation of 
dust abundances and composition with metal enrichment in galaxies spanning a wide range 
of metallicities. Since the metallicity of galaxies evolves monotonically with time, galaxies 
with different metallicities provide snapshots of the evolutionary history of galaxies. 

An exciting result provided by ISO spectral observations of nearby galaxies was the 
discovery of a striking correlation between the strength of their mid-IR aromatic features 
and their metallicity (Madden et al. 2006). Low-metallicity galaxies exhibited very weak 
or no aromatic features. Observations of the 8-to-24 /am bands flux ratio obtained with 
Spitzer/ IRAC and Spitzer/ MIPS instruments showed a correlation of this flux ratio with the 
galaxies’ oxygen abundance (Engelbracht et al. 2005). Since the Spitzer /IRAC^m. band is 
supposed to trace the strength of the aromatic features, and the Spitzer/ MIPS 24 ^ m that of 
the continuum emission from the hot non-aromatic dust component, this correlation seemed 
to confirm the trends discovered by ISO. This correlation was independently confirmed by 
Wu et al. (2006) and O’Halloran et al. (2006) from spectral observations of low-metallicity 
blue compact dwarf galaxies with the Spitzer/ IRS instrument. 

The aromatic features are most commonly attributed to the vibrational modes of Poly- 
cyclic Aromatic Hydrocarbons (PAHs; Leger & Puget 1984; Allamandola et al. 1985), which 
are large planar molecules made of 50 to 1000 carbon atoms. Their ubiquity makes them 
an important component of dust models (Desert et al. 1990; Dwek et al. 1997; Draine & 
Li 2001; Li & Draine 2001; Zubko et al. 2004), locking up ~ 15 — 20% of the total amount 
of interstellar carbon (Zubko et al. 2004, with solar abundance constraints). They mostly 
reside in galactic photodissociation regions (PDRs), where they play an important role in 
the heating of the gas by providing photo-electrons (e.g. Tielens & Hollenbach 1985; Bakes & 
Tielens 1994), and in interstellar chemistry by providing surfaces for chemical reactions. Be- 
cause of their small sizes, PAHs are stochastically heated by the interstellar radiation fields. 
The relative strength of some of the aromatic features depends on their ionization state, 
and varies therefore significantly with the physical condition of the environment (Hony et al. 
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2001; Vermeij et al. 2002; Galliano et al. 2007). Understanding the evolution of PAHs and 
their relation to the global properties of galaxies is therefore extremely important because 
of the complex mutual influences between PAHs, and their ambient radiative and gaseous 
surroundings. 

Several explanations have been offered for the correlation of the intensity of the PAH 
features with metallicity. The first suggests that the trend reflects an increase in the de- 
struction efficiency of PAHs in low metallicity environments Galliano et al. (2003, 2005); 
Madden et al. (2006). Low-metallicity environments are bathed with harder photons than 
our Galaxy, due to the higher effective temperature of their stars, and their young age. 
Moreover, the paucity of dust allows this hard radiation field to penetrate deeper into the 
ISM, compared to high-metallicity systems, selectively destroying the PAH molecules by 
photoevaporation or photodissociation. This explanation is consistent with the mod- 
els of population synthesis and dust evolution of Dwek et al. (2000). Assuming 
that PAHs Eire efficiently destroyed by UV photons in Hil regions, their models 
showed an evolutionary trend of PAH features with time, as the relative contri- 
bution of ionizing OB stars to the galaxy’s SED decreases with time. A second 
explanation has been proposed by O’Halloran et-al. (2006), who suggested that PAHs are 
destroyed by the numerous shocks observed in low metallicity systems. To support their 
proposition, they showed an anti-correlation between the PAH-to-continuum ratio and the 
[Fe n]25.99^m/[Ne Ii]i2.8i/«n line ratio, the latter being supposedly a shock tracer. The problem 
with this explanation is that there is no observational evidence that PAHs are selectively de- 
stroyed in shocks. On the contrary, Reach et al. (2002) showed that, in the shocked medium 
of 3C391, both the PAH features and the underlying continuum disappear. 

All previous explanations attribute the paucity of PAHs to destructive pro- 
cesses that are more efficient in the early stages of galaxy evolution. In contrast, 
Dwek (2005) suggested that the observed correlation reflects an evolutionary trend of the 
sources of interstellar PAHs with metallichy. PAHs and carbon dust are mostly produced in 
asymptotic giant branch (AGB) stars which, unlike massive stars, recycle their ejecta into 
the ISM after a significantly longer time of main sequence evolution. The observed corre- 
lation of PAH line intensities with metallicity is therefore a trend of PAH abundance with 
galactic age, reflecting the delayed injection of PAHs and carbon dust into the ISM by AGB 
stars in their final, post- AGB, phase of their evolution. 

Such distinct evolutionary trends of SN and AGB produced dust with time 
were predicted by Dwek (1998) and more recently by Morgan &; Edmunds (2003). 
In particular, the latter show the evolutionary trend of AGB dust with time. This 
trend can be translated as a trend with metallicity if galaxies approximately share 
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a global “cosmic” star formation history. 

The goal of this paper is to present a detailed evolutionary model to examine whether the 
observed trend of PAH line intensity with metallicity reflects an evolutionary trend of PAH 
abundance with metallicity. This requires the determination of PAH abundance from the 
strength of their aromatic features in the galaxies for which this trend has been observed, and 
the use of a chemical evolution model to follow the change of PAH abundance with galactic 
metallicity (Dwek 1998). The paper is organized as follows. We first present in Sect. 2 the 
sample of nearby galaxies that were considered in our analysis. In Sect. 3, we describe the 
method we used to determine the intrinsic stellar radiation field, and the abundance of PAHs 
as well as the larger grains in these galaxies. In Sect. 4, we briefly describe the chemical 
evolution model used in the calculations,, and compare its results to the abundances derived 
from our SED modeling. The results of the paper are briefly summarized in Sect. 5. 

Throughout this paper we will refer to the solar abundances by Grevesse & Sauval 
(1998), the oxygen number abundance being 12 + log(Q/H)© = 8.83, the Helium and heavy 
elements to gas mass ratios Y© = 0.248 and Z© = 0.017, respectively. Besides, we assume 
that the helium abundance is independent of the metallicity. 


2. THE SAMPLE OF NEARBY GALAXIES 
2.1. Source Selection 

In order to properly characterize the PAH emission, we considered galaxies whose mid-IR 
spectrum has been observed, either with one of the spectrographs onboard the ISO satellite, 
or with the Spitzer/ IRS. We combined ISO samples of starbursts and AGNs (Laurent et al. 
2000), spirals (Roussel et al. 2001), ellipticals (Xilouris et al. 2004), dwarf galaxies (Madden 
et al. 2006), and the low- metallicity sources of the Spitzer sample presented by Engelbracht 
et al. (2005) which were observed by Spitzer/ IRS. Among the ellipticals, only NGC 1399 has 
a mid-IR spectrum. Consequently, our sample includes various types of galaxies and covers 
a wide range of metallicities and star formation activity (Table 1). 

The modeling that will be presented in §3 requires the assembly of data covering the 
stellar as well as the dust emission components, for each galaxy. First, we need most of 
the U, B, V, R, I, J, H, K fluxes, to constrain the stellar spectrum shape. Photometric 
J, H and K bands are available for almost all of our sources, thanks to the 2MASS survey 
(Jarrett. et al. 2003). We rejected the galaxies Arp 118, Arp 236, Arp 299 and NGC 4038, 
for which no B and V bands were reported. Second, the far-IR SED is used to constrain 
the intensity of the interstellar radiation field (ISRF). Hence, we rejected galaxies which 
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have not been detected by IRAS or MIPS, such as HS 0822+3542, Tol 1214-277 and Tol 65. 
In addition, our mass estimates are normalised to the Hi mass. We therefore rejected the 
galaxy IRAS 23128-5919, for which no Hi observation was reported. Finally, we removed 
M 31 from our sample, since its angular diameter is too large to build a consistent observed 
total SED. 

The global properties of the selected sources are presented in Table 1. If relevant, the 
distances were homogenised to Hq = 71 kms' 1 Mpc -1 . The masses have been scaled to the 
adopted distance. 





Table 1. Select properties of the sample. 


Same 

R.A. 

(J2000) 

Dec. 

(J2000) 

Mid-IR 

spectrograph 

Distance 

(Mpc) 

12 + log(0/H) 
[ref.] 

M H i 

(10 s Mq) 

[ref.] 

A/h 2 

(10 8 Mq) 

[ref.] 

Notes 

> 11 

00 ,l 36 m 52!5 

— 33°33 / 19 ,/ 

Spitzer/ IRS 

92 

7.9 

[2] 

< 1(?) 

[2] 

<1 

[2] 

Pec H II 

3 253 

00 h 47 m 32*9 

—25° 17' 18" 

ISO/CAM 

3.3 

9.0 

[3] 

18 

[4] 

17 

[5] 

Hll 

2 520 

01 ,l 24 m 34!9 

+03°47'31" 

ISO /CAM 

27 

(?) 


35 

[6] 

35 

[7] 

Pec H II 

J 613 

01 fc 34 m 17?5 

— 29°24 , 58" 

ISO /CAM 

19 

9.2 

[8] 

37 

[9] 



Sy 

^ 891 

02 ,l 22 m 33®4 

+42°20'57" 

ISO/ CAM 

9.6 

8.9 

[10] 

76 

[9] 

49 

[9] 

Edge-on 

3 1068 

02 h 42 m 40?6 

-00°00'47" 

ISO/CAM 

15 

9.0 

[11] 

22 

[12] 

74 

[13] 

Sy 

J 1097 

02*46 m 19?l 

— 30°16 / 28" 

ISO/C AM 

12 

9.0 

[14] 

51 

[15] 

> 7.3 

[16] 

Sy 

3 1140 

02 h 54 m 33?5 

— 10°01'44" 

ISO/ CAM 

25 

8.0 

[17] 

52 

[18] 



Irr Hll 

3 1365 

03 h 33 m 35?6 

— 36°08'23" 

ISO/ CAM 

19 

9.1 

[19] 

130 

[20] 

170 

[21] 

Sy 

0335-052 

03 ,1 37 7n 44*!0 

-05°02'38" 

Spitzer/ IRS 

54 

7.3 

[22] 

9.9 

[23] 



BCD 

3 1399 

03 /l 38 ,n 29n 

-35°27'03" 

ISO /CAM 

21 

(?) 


< 2 

[9] 



cD 

42 

03 ,l 46 m 49!7 

+68°05'45" 

ISO/CAM 

3.8 

8.9 

[24] 

130 

[9] 

50 

[9] 

Hll 

3 1569 

04 ,l 30 TO 49?l 

+64°50'52" 

ISO/CAM 

2.2 

8.2 

[25] 

1.3 

[26] 

0.50 

[27] 

Irr H II 

3 1808 

05 /l 07 TO 42!3 

— 37°30'47" 

ISO/ CAM 

11 

9.1 

[28] 

18 

[29] 

20 

[30] 

Sy 

tv 40 

05 ,l 55 m 42*7 

+03°23'29" 

ISO/CAM 

10 

8.1 

[31] 

4.2 

[9] 

0.23 

[9] 

BCD 

!-10 

08 fc 36 m lK2 

-26°24'34" 

Spitzer/ IRS 

8.7 

8.9 

[25] 

3.1 

[32] 

1.6 

[32] 

Irr Hu 

r 18 

09 /l 34 m 02®0 

+55°14 , 28" 

Spitzer /IRS 

13 

7.2 

[22] 

1.2 

[33] 



BCD 

2 

09- l 55 m 51®8 

+69°40'46" 

ISO/CAM 

3.6 

9.0 

[34] 

9.0 

[35] 

13 

[36] 

Irr H II 

3 3256 

10 fc 27 w 5Hl 

-43° 54' 17" 

ISO/CAM 

37 

8.9 

[37] 

62 

[38] 

300 

[39] 

Pec Hll 

33 

10 ft 32 m 31?9 

+54°24'04" 

Spitzer/ IRS 

20 

8.4 

[37] 

4.3 

[40] 

0.70 

[41] 

Irr H II 

153 

10 /l 49 m 05*0 

+52°20'08" 

Spitzer /IRS 

37 

7.8 

[42] 

6.7 

[43] 



BCD 

Zw 403 

ll ,l 27 m 59?9 

+78°59'39" 

Spitzer/ IRS 

4.5 

7.7 

[22] 

0.69 

[40] 



BCD 

448 

ll h 42 m 12!4 

+00°20'03" 

Spitzer /IRS 

70 

8.0 

[44] 

47 

[45] 

24 

[45] 

Pec H Ii 

3 4945 

13 /l 05 m 26®2 

-49°28T5" 

ISO/ CAM 

3.9 

(?) 


45 

[46] 

17 

[47] 

Edge-on S? 

taurus A 

13 h 25”*28?0 

-43°01'06" 

ISO/CAM 

3.8 

9 

[48] 

11 

[49] 

2.0 

.[50] 

Elliptical 9. 

1 

13 ,t 29 m 52?7 

+47°11'43" 

ISO/CAM 

8.4 

8.7 

[51] 

50 

[52] 

56 

[53] 


3 

13 h 37 m 0ffi7 

— 29°51'58" 

ISO/C AM 

4.5 

9.2 

[54] 

51 

[55] 

49 

[56] 

H II 



Table 1 — Continued 


Name 

R.A. 

(J2000) 


Mid-IR 

spectrograph 

Distance 

(Mpc) 

12 + log(0/H) 
[ref.] 

Mr i 
(10 8 M e ) 

[ref.] 

M H2 

(10 s Me) 

[ref.] 

Notes 

189 

14 /l 01 m 21!5 

— 33°03'50" 

Spitzer/ IRS 

15 

8.0 

[57] 

13 

[58] 

. . . 


Hn 

rcinus 

14 ,l 13 m 09?6 

— 65°20'21" 

/50/CAM 

4.0 

(?) 


55 

[46] 

11 

[59] 

Sy 

3C 5253 

is^^-ss?? 

— 31°38'29" 

/50/SWS 

3.3 

8.2 

[17] 

0.91 

[60] 

< 0.1 

[61] 

Irr Hll 

p 220 

15 /l 34 rra 57^2 

+23°30'11" 

/50/CAM 

73 

(?) 


280 

[38] 

160 

[62] 

ULIRG 

JC 6240 

16 h 52 m 58®8 

+02°24'06" 

/50/CAM 

98 

(?) 


88 

[9] 

310 

[9] 

LIRG 

1C 6946 

20 h 34 m 51?2 

+60°09'17" 

/50/CAM 

5.5 

9.1 

[63] 

56 

[64] 

33 

[65] 

Hu 

■k 930 

23' l 31 m 58?3 

+28°56'50" 

Spitzer /IRS 

73 

8.1 

[44] 

28 

[66] 



Hn 

1C 7714 

23 /l 36 m 14®l 

+02°09'19" 

Spitzer/ IRS 

37 

8.5 

[67] 

17 

[68] 

22 

[68] 

Pec H II 


3ferences. — [2] Bergvall at ah (2000); [3] Zaritsky et al. (1994); [4] Boomsma et al. (2005); [5] Houghton et al. (1997); [6] Bernloehr 
3); [7] Yun & Hibbard (2001); [8] Alloin et al (1979); [9] Bettoni et al. (2003); [10] Otte et al. (2001); [11] Dutil & Roy (1999); 
Staveley-Smith & Davies (1987); [13] Heifer et al. (2003); [14] Storchi-Bergmann et al. (1995); [15] Ondrechen & van der.Hulst (1989); 
Gerin et al. (1988); [17] Heckman et al. (1998); [18] Hunter et al. (1994); [19] Roy & Walsh (1997); [20] Jorsater & van Moorsel (1995); 
Sandqvist et al. (1995); [22] Izotov et al. (1999); [23] Thuan et al. (1999a); [24] Pilyugin et al. (2004); [25] Kobulnicky & Skillman 
7); [26] Stil & Israel (2002); [27] Israel (1997); [28] Ravindranath & Prabhu (2001); [29] Dahlem et al. (2001); [30] Dahlem et al. 
0); [31] Perez-Montero <fc Diaz (2003); [32] Sauvage et al. (1997); [33] van Zee et al. (1998); [34] Boselli et al. (2002); [35] Appleton 
. (1981); [36] Walter et al. (2002); [37] Mas-Hesse & Kunth (1999); [38] Casasola. et al. (2004); [39] Sargent et al. (1989); [40] Thuan 
. (2004); [41] Israel (2005); [42] Kunth & Joubert (1985); [43] Thuan &: Martin (1981); [44] Izotov & Thuan (1998); [45] Sage et al. 
2); [46] Huchtmeier & Richter (1988); [47] Dahlem et al. (1993); [48] Schaerer et al. (2000); [49] Richter et al. (1994); [50] Wild 
. (1997); [51] Bresolin et al. (2004); [52] Dean & Davies (1975); [53] Heifer et al. (2003); [54] Webster & Smith (1983): [55] Tilanus 
ien (1993); [56] Lundgren et al. (2004); [57] Durret et al. (1985); [58] Paturel et al. (2003); [59] Elmouttie et al. (1998); [60] R,eif 
. (1982); [61] Meier et al. (2002); [62] Sanders et al. (1991); [63] Kobulnicky et al. (1999); [64] Carignan et al. (1990): [65] Tacconi 
>ung (1986); [66] Hopkins et al (2002): [67] Gonzalez-Delgado et al. (1995)r[68] Struck & Smith (2003). 


te. - The sources are ordered according to their right ascension. The symbol identifies uncertain values. Entries for the 




allicity for which no data exist are marked by • • * , and assumed to be solar-to-supersolar. The H2 masses are not used to calculate 
dust-to-gas mass ratio, because the IR emission was assumed to originate entirely from the Hi gas. 
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2.2. ISO/ CAM Data Reduction 

Most of the sources in Table 1 were observed with ISO / CAM (Cesarsky et al. 1996a) 
on board the ISO satellite (Kessler et al. 1996). These spectra were used by Madden et al. 
(2006), and we refer to this work for a detailed description of the data reduction. The CVF 
performed spectral imaging using a 32 x 32 detector array, with a sampling of 3 // pixel~ 1 or 
O^pixel -1 in our cases, from A = 5 pm to 16.5 pm with one pointing of two CVFs, from 
A = 5 to 9.5 pm and from A = 9.0 to 16.5 pm. The spectral resolution goes from A/AA = 35 
to 51 across the full spectra. 

For the data treatment, we used the CAM Interactive Reduction (CIR, version AUG01; 
Chanial 2003). The subtraction of the dark currents was performed using the Biviano et al. 
(1998) model which predicts the time evolution for each row of the detector, taking into 
account drifts along each orbit and each revolution. We masked the glitches using multi- 
resolution median filtering (Starck et al. 1999) on each block of data after slicing the cube. 
Additional deglitching was performed manually, examining the temporal cut for each pixel. 
We corrected the systematic memory effects using the Fouks-Schubert method (Coulais & 
Abergel 2000). We computed a hybrid flat-field image placing a mask on the source and 
computing a flat field outside this mask from the median of the temporal cut for each pixel. 
For the pixels which were on-source, the flat-field response was set to the corresponding 
calibration fiat-field. The conversion from Analog Digital Units to mJy/pixel was performed 
using the standard in-flight calibration data base. To remove the sky contribution, the 
sources smaller than the array were masked and, for a given wavelength, the median of the 
pixels which are off-source were subtracted from each pixel. For the more extended sources, 
we subtracted an independently observed zodiacal spectrum. The intensity of this spectrum 
was a free parameter varied in order to match the properly sky subtracted fluxes in the LW2 
(6.7 pm) and LW3 (14.3 pm) broadbands. The final product is a 3D spectral-image of each 
galaxy. We integrated the spectrum into an aperture encompassing the entire galaxy, and 
to obtain the global SED of the galaxy. When the angular size of the source was larger than 
the one of the array, we scaled-up the spectrum to match the IRAS^m broadband flux. 


2.3. Spitzer/ IRS Spectrum Extraction 

Several of the low-metallicity sources in Table 1 were not observed by ISO/ CAM. 
we therefore complemented our database with public-ally released mid-IR spectra from the 
Spitzer/ IRS spectrometer on board the Spitzer Space Telescope (Houck et al. 2004b; Werner 
et al. 2004). Among these galaxies, the spectra of SBS 0335-052, NGC 7714, VII Zw 403, 
Haro 11 and I Zw 18 are described in details in Houck et al. (2004a), Brandi et al. (2004) and 
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Wu et al. (2006). We considered only low-resolution data, taken with the SL (Short-Low) 
module, from A = 5.2 pm to 14.5 pm, and the LL (Long-Low) module, from A = 14.0 /im 
to 38.0 pm, both with a spectral resolution of A/AA ~ 64 — 128. 

. We retrieved the Basic Calibrated Data (BCD) that have been preprocessed by the 
Spitzer Science Center (SSC) data reduction pipeline, and converted to flux density, corrected 
for stray light and flatfielded. The extraction of the spectra from the 2D space/ wavelength 
images was performed with the Spectral Modeling, Analysis and Reduction Tool (SMART, 
version 5.5.6; Higdon et al. 2004). We first inspected the BCD images and identified the 
hot pixels which have not been masked out by the SSC. We replaced them by the median 
of their neighbors. For each module, an off-source position is observed. We subtracted this 
spectrum from the on-source one, in order to remove the sky emission. The extraction of 
the ID spectrum was performed inside a column whose width vary with the wavelength. 
We have excluded the bonus order. Then, the various frames, for each nod position, were 
coadded. Since the long wavelength end of the SL module and the short wavelength end of 
the LL one were not systematically overlapping, we finally scaled the SL module, in order to 
obtain a continuous spectrum. This scaling factor can be as large as 50%. At the time when 
this publication is written, the Spitzer / IRS data handbook recommends not to derive the 
signal-to-noise ratio from the uncertainties generated by the pipeline. Instead of that, we 
adopted the recommended systematic error of 20%. However, the SL module of the noisiest 
spectra (VII Zw 403, Mrk 153, Mrk 930, I Zw 18, Tol 89) exhibits fluctuations larger than 
this value. To take into account these statistical variations, we smoothed these spectra into 
a AA ~ 0.2 pm window (4 points), and took the standard deviation inside this window as 
the error. Similarly to what we did with ISO / CAM spectra (§2.2), in order to compensate 
the fact that we may be overlooking some extended emission, we scaled-up the Spitzer/IRS 
SL and LL spectra to match the IRAS^m and IRAS^m- In the case of II Zw 40, the 
Spitzer / IRS spectrum does not exhibit the PAH features that Madden et al. (2006) detected 
in the extended emission, thus we prefer to use the ISO /CAM data for this galaxy. For 
Tol 89, the Spitzer /lRAC^ m broadband flux is higher than the integrated spectrum into 
the same band. This is probably due to the fact that the slit of Spitzer / IRS measures only 
the nucleus emission, which likely has a steeper continuum and weaker features than the 
extended emission contributing to the total broadband. Therefore, we will consider the mass 
of PAHs derived for this galaxy to be a lower limit. 
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2.4. The Multiwavelength SEDs 

The UV-to-FIR SED of each galaxy was built using data from different catalogs. We 
extensively used the HYPERLEDA and NED databases. Most of the SEDs include IRAS 
broadband fluxes at 12, 25, 60 and 100 pm (Moshir & et al. 1990; Rice et al. 1988), and 
J, H, K 2MASS fluxes (Jarrett et al. 2003). For some of the galaxies, we completed the IR 
SED with the broadband Spitzer data reported by Dale et al. (2005). We kept only data 
encompassing the entire galaxy. The optical data were corrected for Galactic extinction, 
using the Schlegel et al. (1998) extinction maps. 

For NGC 5253, we used the ISO/ SWS observations of Crowther et al. (1999) to model 
the spectrum of this galaxy. 


3. DETERMINATION OF THE GLOBAL MASS OF PAHS 
3.1. Motivations and Approach 

The 8/24 pm flux ratio provides information on the strength of the PAH features in 
the galaxy, although it is relatively limited. Indeed, the Spitzer /IR AC flux is dominated 

by the 8.6 pm PAH feature, but can be dominated by the continuum, in the case of low- 
metallicity galaxies (Madden et al. 2006), or by the silicate extinction feature at 9.7 pm, 
in the case of deeply embedded sources (e.g. Spoon et al. 2004; Gallais et al. 2004). When 
there are no PAHs and the 8 pm is dominated by the continuum, the 8/24 pm ratio reflects 
the color of the mid-IR emitting silicates and graphites. The SpitzerfMIPS^m flux is 
dominated by the continuum emission from, the silicate and carbon grains. However, this 
part of the SED is very sensitive to the abundance of small grains (e.g. Galliano et al. 2003, 
2005), and can originate additionally from large grains close to strong sources of radiation 
(e.g. Plante & Sauvage 2002; Vanzi & Sauvage 2004). Fig. 1 shows the ratio of the two 
bands, for the sources presented in Table 1. The error bars on the ratio, in Fig. 1, come 
from the propagation of the observational errors. The error bars on the metallicity were not 
systematically given by the authors who published them. An error of 0.1 dex in the O/H 
number abundance accounts for the typical dispersion between independent measurements. 
There were 6 sources, in Table 1, for which no metallicity measurements were reported. 
However, all of them have the morphology, optical colors, IR emission and H 2 content of 
solar or supersolar systems, with NGC 1399 probably having the largest metallicity. We 
assign an arbitrary Z — (2 ± 1) Z 0 to these galaxies, and consider them to be uncertain. 

The optical properties of the PAHs have not the same spectral shape than the ones 
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10 -5 10' 4 10" 3 

0/H 


Fig. 1. — Mid-IR color as a function of the oxygen number abundance, for the sam- 
ple in Table 1. Our observed SEDs have been integrated over the Spitzer/IRAC Sfim and 
Spitzer/ M.W& 2 A fim. bandpasses, in order to produce this plot. The two open circles are 
the (U)LIRGs of our sample. The grey stripe is the ±lcr linear correlation between 
log(0/H) and log (F„(8 pm)/ F u (2 A pm)). The hatched stripe shows the range of the ra- 
tio (0.06 < i\,(8 pm)/F u ( 24 pm) < 0.13), when it measures the color of the silicate and 
graphite continuum due to the weakness of the PAH features. 
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of solid-state grains (carbonaceous or silicates) believed to be the carriers of the mid-IR 
continuum. Fig. 2 shows the wavelength-dependence of the mass absorption coefficient, « a b S , 
of various types of grains. It shows that the value of «; a b s for PAHs drops by 4 orders of 
magnitude between 0.1 and 1 fim , where the stars emit most of their energy, while the « a bs 
of graphite and silicate dust drop by less than 2 orders of magnitude. However, most of the 
energy is absorbed at shorter wavelengths. Thus, the PAHs are more sensitive to very young 
stellar populations than the grains responsible for the 24 [xm continuum (Fig. 3), and a 
correct determination of their excitation rate should ideally take into account this property. 
Beyond 1 /xm, the opacity of analog PAHs measured in laboratory (Mattioda et al. 2005) 
can be higher than the values stated by Li & Draine (2001). 



Wavelength [/xm] 

Fig. 2. The wavelength-dependence of the mass absorption coefficient of the PAHs (Li & 

Draine 2001), graphites (Laor & Draine 1993) and silicates (Weingartner & Draine 2001), 
integrated over the Zubko et al. (2004) size distribution, for the bare grain model with solar 
abundance constraints (BARE-GR-S). 
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Wavelength [/i,m] 

Fig. 3. Sensitivity of the energy absorbed by the PAHs to the age of the stellar populations. 

The solid lines are two instantaneous bursts modeled with PEGASE, with a Salpeter IMF. The 
black spectrum is taken immediately after the burst, while the grey one has evolved during 
1 Gyr. The dashed lines are the part of these spectra absorbed by the PAHs, with the 
mass absorption coefficient of Fig. 2, and a thin slab geometry. The hatched areas are the 
energies carried by non-ionizing photons absorbed by the PAHs, for each stellar spectrum. 
This energy is 10 times higher for the young burst, than for the old one. 
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The determination of the mass of PAHs, from an observed SED, requires a more sophis- 
ticated approach than a simple modified black body fit, since they are transiently heated, 
hence only a fraction of them is excited at a given time. Consequently, a proper determi- 
nation of the PAH content must also take into account the characteristics of the intrinsic 
radiation field, in addition to their size distribution. 


3.2. Determination of the Global Interstellar Radiation Field 

In order to study the variations of the PAH content with the metallicity, we considered 
several galaxies at different metallicities, some of these galaxies being scarcely observed. 
Thus, we developed a simple self-consistent model, fully constrained by the data. The major 
assumption of our approach is that the PAHs are bathed by an homogeneous radiation field 
throughout the galaxy. More precisely, we assume that PAHs are not present in H II regions, 
where they are destroyed (e.g. Cesarsky et al. 1996b; Peeters et al. 2002; Martm-Hernandez 
et al. 2002), and that they are not present in molecular clouds, where they are probably 
accreted onto larger grains. In star forming regions, the gas density to radiation 
field intensity is generally such that the H 1 /H 2 transition lies very close to the 
surface and that the gas in the UV penetrated cloud layer is partially molecular. 
However, most of the molecular mass is devoid of PAH emission. We therefore 
assume that all the PAHs are associated with the H I gas, and that they do not receive any 
ionizing photons (A < 0.0912 /zm), which are absorbed into the Hil regions. Finally, we 
assume that the grains emitting most of the far-IR energy are bathed by the same radiation 
field as the PAHs. This latter assumption is corroborated in part by the observed correlation 
between the intensity of the PAH 7.7 pm feature and the 850 fj,m flux in galaxies (Haas et al. 
2002). However, this correlation may not be valid at spatial resolutions high enough to resolve 
the clumps. 

The determination of the intensity of the ISRF is done in the following steps. 

Determination of the intrinsic stellar luminosity. We use the stellar evolution and 
population synthesis model PEGASE (Fioc & Rocca-Volmerange 1997), in order to gen- 
erate a proper stellar spectrum, for each galaxy. Its luminosity, L*, is the linear 
combination of two instantaneous bursts, of different ages. The stellar monochromatic 
luminosity is: 

= Afyoung X T^tyoung) + ^o\d X L\,(t, 0 ]d) , (1) 

M* j d and M* oung being the masses of each components, and t 0 w and t you ng their ages. 
These are the free parameters. L*(£) is the luminosity per unit mass, t Myr after 



-16- 


an instantaneous burst, assuming a Salpeter initial mass function. These two bursts 
represent the old underlying stellar population and the most recently formed one, 
respectively. This modeling does not take into account a possible contribution from an 
AGN, since in most of our galaxies, the SED is dominated by the stellar input. 


Estimation of the escaping stellar light fraction. The escaping stellar radiation is: 


^UV-opt. = X 'P esc ('T, w), 


( 2 ) 


where r(v) is the optical depth of the medium and u(u) the albedo of the dust. The 
escaping fraction, 'P esc , depends on the geometry of the medium. We consider consec- 
utively two extreme geometries: a sphere where stars and dust are uniformly mixed 
(Varosi Sc Dwek 1999): 


■P e8 c(r, to) 


p e(r) 

1 — w[l — P j(t)] ’ 


(3) 


with 




(4) 


and a thin shell geometry with central illumination, for which: 


VesciTiLo) = e 


(5) 


Another unknown is the shape of the extinction law, r(v), and its corresponding albedo, 
to{v). In order to avoid strong hypothesis on the dust properties, we consider consec- 
utively Galactic, LMC and SMC extinction curves. Similarly to what Pei (1992) did, 
we fit these three curves with a linear combination of graphites and silicates, the PAHs 
representing only a small fraction of the mass are neglected at this point. We adopt 
the BARE-GR-S grain size distribution (bare grains with solar depletion constrains 
on the amount of metals locked into the dust; Zubko et al. 2004). The opacity of the 
grains is: 

K'dust^) = Cgra X Kg ra (t / ) + Gil X ^sil^)) ( 6 ) 

where K gl . a (i^) and k s h(u) are plotted at Fig. 2, and their derived relative mass fraction, 
Cgra and Giii are given in Table 2 . For an homogeneous medium, with a dust mass 
density of pduso 

t{v) = Kdust (") x / At ust ds (7) 

Jo 

= *dust (v) X A, (8) 


where L is the line of sight and A is the dust mass column density. Thus, for a given 
extinction type and geometry, the only free parameter governing V esc (r, u>) is A. The 
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Table 2: Mass fractions of silicate and graphite, for different extinction curves. 



Galaxy 

LMC 

SMC 

<r 

0.36 

0.14 

0.02 

cf 

0.64 

0.86 

0.98 


escaping luminosity, Z^ v ' opt ', is constrained' by the observed UV-optical broadband 
fluxes. 


Derivation of the intrinsic radiation density. The intrinsic stellar monochromatic lu- 
minosity gives us the global shape of theTSRF. The radiation energy density, U v , is 
defined as follows: 


L*(z/) x [1 - P es c(t,u>)]. = c JjJ^pdust(r)K dust (u) x U v {f, u) d 3 r. 


( 9 ) 


The left-hand side is the global absorbed energy, and the right-hand side is the inte- 
gral over the volume of the locally absorbed energy. In the case of an homogeneous 
distribution, it gives: 


U v {y) 


' -T i x L »< 
cM dust K dust {U) 


(10) 


where Md US t is the total - non PAH - dust mass. This radiation density is constrained 
by the far-IR emission. Indeed, this part of the SED is not sensitive to the size 
distribution, since it is dominated by large grains at equilibrium with the radiation 
field, still small enough to be in the Rayleigh regime. Thus, we adopt the BARE-GR-S 
of the Zubko et al. (2004) model and fit the far-IR SED (usually the IRASeo^ m and 
IRAS m ^ m bands) with: 


L™ d “» = Aftm d... x (Cgr,4‘ r » + Chi 4*"M) , (11) 


where Lt sra< ' h ' 1) is the graphite (silicate) monochromatic luminosity per unit mass of 
dust. Since, the grains emitting in the far-IR are the largest, and they contain the 
essential of the mass, we assume that: 


A/dust — AfpiR dust • (12) 

We will discuss, in §4, possible deviations to this relation. 


Fig. 4 shows an example of a SED fitted in order to determine the radiation den- 
sity. The B, V, H, K bands, at least, are required to constrain the burst parameters, 
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Myoungytyovmg, M 0 id, told- To constrain the parameters of the ISRF, A and Md US t, we ad- 
ditionally need at least the 60 and 100 fj,m IRAS fluxes. The excess between ~ 10 and 
~ 40 fj,m can be attributed either to an overabundance of small grains or to hotter regions, 
while a possible excess in the submillimeter could be attributed to colder shielded regions. 
We do not require energy conservation between the energy absorbed and reradiated by the 
dust, since we fit the SED only partially. 



Wavelength [/i.m] 


Fig. 4. — Determination of the intrinsic stellar radiation density. The symbols and the mid- 
IR spectrum are the observations. The mid-IR regime is purposely not fitted at this step. 
The open circles are the data which have not been used to constrain the ISRF. 


3.3. Estimation of the PAH Mass 

The radiation density in Eq. (10) is used to excite two components of neutral and 
ionised PAHs, with the BARE-GR-S size distribution of Zubko et, al. (2004) and the Li & 
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Draine (2001) optical properties, taking into account the stochastic heating. The luminosity 
reemitted by the PAHs is: 

P AS (v) = M pah x ([1 - /+] £» PAH V ) + ULt rAH *(v)) , (13) 

where Mpah is the mass of PAHs, /+, the mass fraction of cationic PAHs, and £,t pAH °( PAH+ ) ) 
the luminosity per unit mass of PAHs. The free parameters controlling this component are 
/+, and Mpah- We incorporate several other physical ingredients in order to fit the detailed 
observed mid-IR spectrum. A very hot grain (VHG) continuum is modeled with the sum of 
two modified black bodies, with the absorption efficiencies of graphite and silicate, Ql hs (a, v). 
and their densities pi. This approximation can not be valid in the case of stochastically 
heated very small grains, however this analytical formulation is sufficient to subtract the 
underlying continuum present in all our spectra. A more sophisticated component would 
be underconstrained. The silicate feature in emission is usually not prominent, except in 
Mrk 153. The luminosity of this VHG continuum is: 


- VMCJ 


= E 


3 irK'"° /QibsK*')' 


a o 


x 


(14) 


for i = { graphite , silicate}. The masses M™ 0 , and the temperatures 7) VHG are the only free 
parameters. The grain radius oo is arbitrarily fixed to 0.03 pm, the results being insensitive 
to its value. 

The mid-IR contains several powerful emission lines. The brightests are usually the 
[Ar n] 6 . 98 ^m> [Ar in]g,gg^ m , [S iv] io. 5 i^im > [Nh and [Ne III] i 5 . 56 ^m • ^Ve therefore add to 

the fit our mid-IR spectra N ion ionic lines of luminosity U° n , with gaussian profiles: 


Ns. 


LT 


V2 


(u) = yL i -^=- 
K W^Av 

i=i v 


x exp 


-2 


V -Vi 

Av 


(15) 


The width, Av, is set by the spectral resolution. The only free parameters of this component 
are the various L,. These lines are not used in this paper, we model them in order to properly 
subtract their contribution to the total spectrum. 

Finally, all these components are attenuated by the escaping probability, V esc . In some 
cases, where the sources are deeply embedded, the optical depth deduced from the UV can 
differ significantly from the one derived from the observed silicate feature absorption at 
9.7 pm. This effect has been encountered in numerous ISO spectra (e.g. Thuan et al. 1999b; 
Gallais et al. 2004). Therefore, we determine another optical depth, r MIR (^), from the fit of 
the mid-IR spectrum. 
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In conclusion, the spectrum is modeled by the following equation: 

M = K(v) x Pes C (r, w) + + 4“(0) x 7V(r MIR , w), (16) 


the stellar continuum being important only for elliptical galaxies. Fig. 5 demonstrates this 
method. From this fit, we derive the total PAH-to-total-gas mass ratio, Z^xa'- 

Zpah = ^ PAH where n = — (17) 

PAH /im h , ^ i-y-z’ v 1 

and Z is the metallicity of the galaxy. The mass fraction of carbon which is locked in PAHs 
is: 

Zq in PAHs — ZpAH • (18) 
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[Aril] [Arlll] [SIV] [Nell] [Nelli] 



Rest frame wavelength [/xm] 


Fig. 5. — Determination of the PAH abundance. The error bars are the observed mid-IR 
spectrum. The dashed line is the sum of the PAH 0 and PAH + spectra. The dashed-dotted 
line is the intrinsic hot grain continuum. The black solid line is the attenuated sum of all 
these components, fitted to the observed spectrum. 
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4. THE PAH ABUNDANCE IN GALAXIES 
4.1. Results of the SED Modeling 

We applied the method described in §3.2 and 3.3, to the sample in Table 1. For each 
galaxy, we performed the SED and mid-IR spectrum fitting for the three dust compositions 
given in Table 2, and for the two geometries: uniform and thin shell. We show the figures 
obtained with Galactic extinction and uniform geometry, in Figs. 6 to 14. The other models 
give fits of the same quality, the main difference being due to the variations in the strength of 
the 9.7 jxm silicate feature in extinction. The dispersion of the results given by the different 
models is an indication of their sensitivity to the model parameters. Consequently, the upper 
and lower limits presented in Fig. 15 and Table 3, are the maximum and minimum value of 
each quantity derived by the six models. In addition, these maximum and minimum values 
take into account the observational errors (e.g. the upper limit on the mass of PAHs is the 
value of the mass given by the model which gives the higher mass, when fitted through the 
higher end of the error bars on the observed spectrum). In most of the cases, the dispersion 
among the different models is larger than the observational error. 

PAH features are not detected in I Zw 18, SBS 0335-052, Mrk 153, NGC 5253 and 
NGC 1399. For these galaxies, we get an upper limit by fitting the short-wavelength part of 
the mid-IR spectrum with the maximum PAH amount allowed by the observational errors. 

I Zw 18: The Hi extends far out of the star forming region. Instead of normalizing the 
dust masses by the value of Mm in Table 1, we consider that the Hi mass associated 
to the star forming region is the one of the- object Hi-A (van Zee et al. 1998), where 
M Hi = 4:4 x 10 7 M©. 

SBS 0335-052: The 65 pm flux is the one reported by Hunt et al. (2005). The fit of this flux 
gives a far-IR dust temperature slightly colder than what we would obtain by fitting 
only the slope of the Spitzer/IRS spectrum. However, this gives a conservative solution, 
since we derive only an upper limit on the mass of PAHs, in this galaxy. The H I halo 
extends also far out of the star forming region. We assume that the star forming region 
has the same size than the optical galaxy, i.e. a radius of 1.6 kpc (6"). We derive the 
corresponding mass of Hi, using the average column density of Am = 7.4 x 10 20 cm -2 
reported by Pustilnik et al. (2001). We find Mm = 4.6 x 10 7 M©. 

VII Zw 403: This Spitzer/IRS spectrum of this galaxy has been studied by Wu et al. (2006). 
However, they did not report any PAH detection. The degradation of the spectral 
resolution that we performed on this spectra (see §2.3) increases the signal-to-noise 
ratio significantly. We report a 4cr detection of the 7.7 pm feature, and marginal 
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detections of the 6.2 and 8.6 pm features. This is the lowest metallicty PAH detection 
to date. 

Haro 11: As quoted by Bergvall et al. (2000), the ratio Mhi/Lb < 0.01 is abnormally low 
in this galaxy. The low value of the H I mass explains the high value of its dust-to- 
gas mass ratio (Fig. 15). On the contrary, its PAH-to-far-IR-dust mass ratio is not 
peculiar. This source, at 92 Mpc, is one of the most distant object in our sample. 

IC 342: This galaxy is located at 10° of Galactic latitude. Hence, it is highly extracted by 
the foreground, which explains the peculiar shape of its optical/near- IR observed SED. 

Tol 89: The mid-IR spectrum of this galaxy samples only its nucleus. That is the reason 
why it underestimates the Spitzer/lKAGs^m flux. To be conservative, we derive the 
lower limit on the PAH mass by fitting the spectrum, and the upper limit, by fitting 
the broad band. 

NGC 1399: This object is a cD galaxy. Its interstellar medium is likely very tenuous and 
the mid-IR emission is dominated by the contribution of evolved stars. It has not been 
represented on the upper panel of Fig. 15, since we only know the upper limits on the 
PAH and Hi gas masses. However, it is shown on the lower panel; the higher open 
circle is the ratio between the upper limit on the PAH mass and the lower limit on the 
dust mass. 

Circinus: This galaxy is located at Galactic latitude below 5°. Like IC 342, its optical 
fluxes are very uncertain. 

NGC 6240 and Arp 220: These two objects are the (U)LIRGs of our sample. In partic- 
ular, the PAH-to-neutral-gas mass ratio of the most luminous, Arp 220, is very low. 

The top panel of Fig. 15 shows the variation of the PAH and FIR-dust to gas mass 
ratios with the metallicity of the interstellar medium. Each individual galaxy can be seen 
as a snapshot of galaxy evolution, at a given time. First, we note that the Galactic values 
of the dust-to-gas mass ratios are in agreement with the one of the other galaxies, around 
the same metallicity. These Galactic values were derived by Zubko et al. (2004), from the 
fit of the emission and extinction of the diffuse interstellar medium, with further constraints 
from the elemental depletion pattern. Thus, it is a very reliable estimation. This comparison 
confirms that our method does not overlook a significant amount of dust, at least around 
the solar metallicity. Second, the trends of Zpah and ^Fin-dust; with the metallicity are not 
identical. Our sample spreads two orders of magnitude in metallicity. We can see that the 
PAH-to-gas mass ratio rises by five orders of magnitudes, while the FIR-dust-to-gas mass 
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ratio, by only three. This differential evolution is illustrated in the lower panel of Fig. 15, 
showing the PAH-to-FIR-dust mass ratio. This figure is the analog of Fig. 1, but instead of 
considering integrated fluxes, it deals with abundances. The PAH-to-FIR-dust mass ratio 
rises by two orders of magnitude, over our sample, while the IRAC 8 ^ m /MIPS 24 Aim band ratio 
varies only by one order of magnitude. 
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Fig. 6. — Fit of the galaxies’ SEDs. These figures are the analogs of Figs. 4 and 5, for 
all the galaxies. We show only the model for a Galactic extinction law and a uniform 
sphere geometry. The left panel shows the fit of the total SED. The normal lines are the 
components of the model, the escaping stellar and dust emission being plotted with solid 
lines, and the intrinsic stellar with dashed lines. The circles with error bars are the broad- 
band observations, and the thick line is the observed micl-IR spectrum. The right panel 
shows the detailed fit of the rnid-IR spectrum. The grey region is the observed spectrum 
(white line) ±la, the grey dotted line is the stellar continuum, the grey dashed line, the sum 
of cationic and neutral PAH emission, the grey solid line, the very hot grain continuum, and 
t he thick line, the attenuated sum of all the components. For clarity, we did not represent 
the ionic lines, although they are part of the final fit. 
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ig. 7. — Fit of the galaxies 5 SEDs. See Fig. 6 for details. 
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Fit of the galaxies 1 SEDs. See Fig. 6 for details. 
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Fig. 11. — Fit. of the galaxies' SEDs. See Fig. 6 for details. 
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ig. 13. — Fit. of the galaxies' SEDs. See Fig. 6 for details. 
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Fit of the galaxies’ SEDs. See Fig. 6 for details. 
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PAH and far-IR dust masses 



Fig. 15. Top panel : dust-to-gas mass ratio for PAHs, 2 Tpah> and the dust that gives rise 

to the far-IR emission, ZpiRdust) as a function of metallicity. Bottom panel : mass ratio of 
PAH-to-far-IR-dust, as a function of metallicity. The circles correspond to galaxies and the 
open stars to the the diffuse Galactic ISM (Zubko et al. 2004). The filled circles are the 
reliable measurements, and the open circles are the ones which are considered uncertain. 
The grey stripes are the ±lcr linear correlation, in logarithmic scale. These figures show the 
two distinct evolutionary trends of Z p A ji and ^FiR-dust with metallicity. 




Table 3. 


-lame 

z 

m pah 

/+ 

■^dust 


(Zq) 

(M e) 


(Me) 

18 

0-023t°;$s 

< 1.3 


200 ±60 

0335-052 

0.030i°;°o! 

<0.5 


740l|®° 

]w 403 

0.072i°;°i| 

0 OQ3 +0,0 ^ 4 
oig 

2-4+ks 


l.2lg;| x io 3 

153 

<60 


(1.8 ±0.5) x 10 4 

11 

n I c>+0.03 
u - lz -0.02 

8.7±|;| x IO 3 

25i®s% 

2-2lg;| x 10 6 

1140 

Q 1 k+ 0.04 

U.IO_ 0 03 

6.8i^ X 10 3 

41+1% 

1.3l°;I x 10 6 
5.6l§;f x 10 6 

L48 

|-J 1 r + 0.04 
U.iO_ 0 03 

4.2±J : | x 10 4 

62lfg% 

9 

Q 1 *+0.04 
u - iO -0.03 

0.8ig | X 10 3 


3.7ll;| x 10 5 

330 

0.17 ±0.04 

2-5±i;| x 10 3 

50l£% 

1.7lg;l x 10 6 

40 

n i q+o.05 

U ’ 18 -0.04 

4.6111 X 10 2 

9 9th% 

(1.9 ±0.6) x 10® 

5253 

0 21 +0 * 06 
u * /A -0.05 

0.23l°-°® 

*' <200 


7.4ll;| x 10 4 

1569 

(2.9 ± 1.6) X 10 2 

3911% 

8.9ll;g x 10 4 

$3 

Q-M±am 

l‘®i L3 X l ° 4 

38ll°% 

5.9lg;l x 10 s 

7714 

oa7±» : \i 

3.6i 2 ;^ x 10 4 

39lf% 

4.4lg’l X 10® 


°- 7 4l8:$ 

(1.1 ±0.6) x 10 6 

44 ± 1% 

2.51H X io 7 
1.6l“;l X 10 7 

2 

i q+ 0.3 
i * u -0.2 

1.41H x 10 6 

5611% 

891 

■i ^+0.3 
A “"— 0.2 

1.5tJ;'“ x 10® 

5211% 

4.9ll;| X 10 7 

10 

1 2 +0 * ;i 
A - Z -0.2 

6 - 5 -3 0 X 103 

4611% 

S-lllis X 10 5 

3256 

•i o-f-0.3 
A * Z -0.2 

l-lioia x 10® 

53l|% 

5.3l 2 ;! x IO 7 

1068 

1 4+ 0 - 4 

4.4±f 4 6 x 10 s 

4811% 

2 -2lo)i x 10 7 

253 

i 4+0>4 
A ' 4 -0.3 

6.4^4 2 x 10 s 

5911% 

3.5ll;J x 10 7 


1 5 +0 * 4 
A '°-0.3 

2.4+?;| x 10 s 

5711% 

4-5lo)i X 10® 

1097 

1 «+°-4 
t'V-O.Z 

X 10 6 

43+1% 

5.611;! x IO 7 

6946 

A * y -0.4 

l-7toJ X 10® 

3611% 

4.811’! X IO 7 

1808 

o n+ 0 - 5 
z * u -0.4 

i.3i°: 7 0 x io 6 

33l 2 % 

2.ill:'l x io 7 

520 

~ 2 ± 1< ? > 

7.9io 6 X 10 5 

8lll% 

2.8lg;l x 10 7 

1399 

~ 2± 1 (? > 

< 1.0 x 10®C ? > 


>4.5 X 10 s ( ? ) 

4945 

~2±1< ? ) 

4.8lf;| x 10® 

73 ±1% 

2 -0— o'l X 10 7 

us 

£8 2 ± 1< ? > 

3.8+| 4 x 10 s 

43 ±1% 

1.9l!;® x 10® 

3240 

a 2± 1 (? ) 

2.1+!;! x 10 6 

401?!% 

1.2 ± 0.3 X 10 s 

20 

~2±1< 7 > 

4.3l 22 i 4 x 10 s 

97ll 9 % 

2 - 7 li:S x io 8 

1365 

ry i +0.5 
Z ‘ A -0.4 

7.4tl l x 10 6 

53llo% 

1.9lg : | x 10 s 

arus A 

ry 1 +0.5 
za -0.4 

6.5^ 2 ' 2 x 10 s 

3011% 

i.ilgii x io 7 

313 

ry rt + 0.6 

Z - Z -0.5 

1.5+ 3 ;2 x 10 6 

38 ± 1% 

3.0l®;l x IO 7 



300 ± 50 

5301!° 

13011! 

160 ± 30 

5.4lg;| x IO 3 

5.811’j x 10 3 

1.0±o! 3 x 10 4 

l-Oloil X 10“ 

1.1 ± 0.2 x 10 3 

1.6 ± 6.7 X 10 3 

137 ±21 

480 ± 90 

360l 4 ” 

7 401318 


T.2±g;| X 10- s 7.2j:|;J x 10~ 4 

].5l®; 4 x 10- 4 S.7±|;f x 10~ 4 

2.1±f;| x 10- 4 9.3 t\°d l x IO" 4 

2.4+^3 4 x 10- 3 1.7ij; 4 x ID" 3 

1.2f{5;g x 10“ 4 3.0±Tt x 10-4 

2.7l} 1 g° x 10- s 5.8lfj x 10~ 4 

5 ' 7 ^ 3 2 8 7 x 10- s 2.7!f;jj x 10- 4 



-36- 


4.2. A Model for the Gas and Dust Evolution 

We have developed a one-zone single-phase chemical evolution model, to follow the 
abundances and composition of the dust and the metallicty as a function of time, in order 
to interpret the results of our SED modeling. In this section, we will give a brief description 
of the physical processes which are incorporated. A more detailed discussion can be found 
in Dwek (1998) and Dwek et al. (2007). Assuming a common star formation history, we will 
use this model to interpret the observed evolutionary trend of galaxies’ SED with metallicity, 
on a global scale. 


4-2.1. Metal Enrichment and Gas Evolution 

In the present paper, we consider a closed box model. We consider the delayed injection 
of material by different stellar progenitors, but we assume that the mixing of the elements 
in the ISM is instantaneous. We adopt a Salpeter initial mass function IMF, <p(m), defined 
as follows: 


4>(m) oc m 2,35 


for 

normalised to 
with 


mi < m < m u 



dm = 1 


f mi = 0.1 M q 

\ m u = 100 Mq 


(19) 


where m is the mass of individual stars. We define the average stellar mass: 

(m) = I m4>(m) dm. (20) 

J mi 

The evolution of the gas mass surface density, E gas (t), with the time t , is: 

= -S SF r (t) + f E sfr ( t - r(m)) <f>{m) dm, (21) 

di Jmt \ m ) 

where Esfr (t) is the mass of star formed per unit time and per unit surface area, r(m), the 
lifetime of a star of mass m, and m e j(m), its mass of gas returned to the ISM. The first 
term of the right hand side of Eq. (21) is the amount of gas removed by the star formation, 
and the second term, is the delayed injection of the gas (H, He and metals), by the various 
progenitors. The ISM metallicity, is defined as: 


ZlS\l(t) = Zaas(t) + Z dust (t), 


( 22 ) 
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Note. — Most of these quantities are defined in §3.2 and 3.3. is the total mass of PAHs (Eq. 13); is the mass 

fraction of cationic PAHs (Eq. 13); M^ ust is the mass of dust dominating the far-IR regime (Eq. 11); is the bolometric 

luminosity emitted by the dust; L star is the intrinsic bolometric luminosity emitted by the stars. The sources are ordered 
according to their metallicity. The symbol W identifies uncertain values. 
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where Z gas (t) is the metal-to-gas mass ratio, and Zdust(i), the dust-to-gas mass ratio that 
will be discussed at §4.2.2. Its evolution is analog to Eq. (21): 

- - -Zisu(t) Esfr (t) + J Ssfr (t - r(m)) 4>(m) dm, (23) 

where Yz{m) is the yield of elements heavier than He, by the stars of mass [m, m + dm]. 


The elemental yields, Yz(m), for the low mass stars (to < 8 M©) are taken from 
Karakas & Lattanzio (2003a, b), and from Woosley & Weaver (1995) for the high mass stars 
(to > 8 M©). Furthermore, a prescription for the star formation rate is required, in order 
to solve these equations. This prescription is given by the Schmidt (1959) law, with the 
coefficients derived by Kennicutt (1998): 


Ssfr OQ 

1 M q yr -1 pc~ 2 


= (2.5 ± 0.7) x 1(T 10 


/ S gas (t) \ 
Vi M©pc- 2 y 


1.40±0.15 


(24) 


and is used to replace EsfrOO i* 1 Eq. (21). Fig. 16 shows the evolution of the total metallicity, 
and of the reduced gas mass ^ gas (£) = E gas (t)/Eg as , where = E gas (0) is the initial gas 
mass surface density. 


4.2.2. Dust Formation and Destruction 

The evolution of the mass surface density, E;(t), of a given dust specie, is the sum of 
three contributions: (i) the rate of dust destruction by star formation, (it) the rate of dust 
condensation in stellar progenitors, (in) the rate of dust destruction in the ISM, by SN blast 
waves: 

—jj - = —Z\ (t) E SF r(£) + f • r^- r - Ssfr(^ — r(m)) &(m) dm ! ^77T> (25) 

J mi {m} ^dustvv 

where Y\{m) is the yield of the considered dust specie by stars of mass [m, m + dm], and 
Tdust (t), the dust lifetime in the ISM. 

The dust yields are derived from the elemental stellar yields described at §4.2.1, following 
Dwek (1998). For low mass stars (to < 8 M©), the dust yields depend on the value of the 
C/O ratio. We assume that the ejecta is microscopically mixed, so that all the excess carbon 
is locked-up in dust, if C>0. If C<0, then we combine all the available Fe, Si, Mg, Ca and 
Ti, with one O atom to produce silicate dust, and titanium oxydes. For high mass stars 
(to > 8 M©), we assume that the ejecta is only maeroscopically mixed, so that both carbon 
and oxygen rich dust can condense. We assume a condensation efficiency of unity for all 
dust species. 




Time [Myr] 


Fig. 16.-— Evolution of the total metallicity (in units of Z 0 ). and of the reduced gas mass 
(Mgas)i for two different initial surface gas mass densities, 
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The lifetime of dust, T^ ust (t) is directly related to the SN II rate (Dwek & Scalo 1980; 
McKee 1986): 


Tdust (t) 


Sgas(^) 


(26) 


-^SN II CO (wism) ’ 

where i?sNn(0 Is the SN II rate per unit area, and (toism) is the average effective mass of 
gas, swept up by a single SN remnant, where the dust is returned back to the gas phase, by 
either thermal sputtering, or grain-grain collision (Jones et al. 1996). The value of (wism) is 
largely unknown. We explore its effects, by varying it from (mjsu) — 0 M 0 (no destruction), 
to (misM) == 300 M 0 (typical destruction). The latter value is typical of our Galaxy. Indeed, 
if ~5x 10 9 Mq is the total mass of gas in our Galaxy, R^ 1 ~ 1/30 yr -1 , its average 
SN rate, and ~ 5 x 10 8 yr, the mean lifetime of an ISM dust particle (Jones 2004), then: 

M Gal 

<">ism> = Wrfcr - 300 Mq. (27) 


oGal rpGal 
^SN ’dust 


Fig. 17 shows the evolution of the total dust content formed by massive stars, and the 
carbon dust formed by AGB stars, for various initial conditions, and destruction efficiencies. 
The SN II dust evolves almost linearly with the metallicity, in absence of destruction, because 
the metal enrichment is dominated by massive stars. In contrast, the AGB carbon dust starts 
rising when the metallicity of the ISM is around 1/20 Zq. This value corresponds to a time 
of ~ 100 Myr (Fig. 16), which the lifetime of the most massive AGB stars. Hence, the carbon 
dust produced by AGB stars is injected into the ISM, with a delay which corresponds to the 
lifetime of the stars. This evolutionary trend was previously noted by Dwek (1998) 
and Morgan &: Edmunds (2003). The change in the slope of the AGB carbon dust, 
around 1 Zq, is simply due to the fact that AGB stars of lifetime longer than ~ 1 Gyr are 
oxygen rich. The dust destruction effects the evolution for Z > 0.1 Zq. 


4.3. PAHs and the Delayed Injection of AGB Carbon Dust 

Prom an observational point of view, PAHs are believed to form in the eireumstellar 
enveloppes of carbon rich AGB stars, and to be subsequently ejected into the ISM through 
stellar winds. However, the direct indications of this formation process, on individual stars 
are very few (Hony et al. 2001; Boersma et al. 2006). The main reason for this lack of 
observational evidence, is probably that the PAHs are not sufficiently excited by their UV 
deficient stellar progenitor. In what follows, we will assume that PAHs are only 
formed in the envelopes of AGB stars. 

Fig. 18 shows the comparison between the dust and PAH abundances derived from 
the observations (Fig. 15), and the ones produced by the evolution model (Fig. 17). The 



dust 
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^gas/^O 


Fig. 17. — Evolution of the dust content with the metallicity of the ISM. We adopt (toism) = 
300 Mq, for the curves with grain destruction. 
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agreement between the PAH-to-gas mass ratio, and the carbon dust production by AGB stars 
is very good. Hence, the peculiar evolution of the PAH content in galaxies can naturally be 
explained by the delayed injection of carbon dust into the ISM, by AGB stars. 

The galaxy He 2-10 is the only outlier. This object is a solar metallicity blue compact 
dwarf galaxy. Its SED is similar to lower metallicity galaxies (Galliano et al. 2005), however 
its metallicity is very uncertain; for example Vacca & Conti (1992) quoted 1/6 Zq. This 
uncertainty may reside in the fact that this galaxy has two cores, resulting from a merger. 
The properties of these nuclei are different in terms of dust absorption (Phillips et al. 1984), 
molecular gas content (Baas et al. 1994), compact source distribution (Cabanac et al. 2005), 
and mid-IR spectrum (Martm-Hernandez et al. 2006). Thus, our global approach may not 
apply to this object. 

Two processes, that were not taken into account in our model, may poten- 
tially have an impact on the observed trend of PAH abundance with metallicity: 
(i) an additional, non-AGB, source of PAHs, and (ii) a radiative destruction 
mechanism, in addition to the mechanical one by SNRs. We argue below, that 
none of these processes will affect the observed trend of PAHs with metallicity. 

In principle, PAHs can form by some interstellar processes, for example, by 
the hydrogenation of small carbon grains. Initially, their abundance will then 
follow the evolutionary trend of SN-condensed carbon dust. However, their 
absence in low metallicity systems suggests that they are efficiently destroyed, 
presumably by shocks or UV radiation, as envisioned in by Madden et al. (2006) 
or O’Halloran et al. (2006). At later times, the rate of PAH production by these 
interstellar processes will follow the evolutionary trend of AGB stars, since they 
will be the major source of carbon dust in the ISM. 

Our model ignored the destruction of PAHs by UV photons capable of sub- 
liming or dissociating these molecules. However, since all hard UV photons are 
produced by massive stars that eventually become SN, the PAH destruction rate 
by UV radiation will follow the same time dependence as their destruction rate 
by SNR. The main difference between the two destruction mechanisms is that 
the PAH destruction efficiency by UV radiatin will be more efficient at earlier 
times, when the dust abundance is low and the ISM more transparent to UV 
radiation. The main effect of the UV radiation will therefore be to destroy any 
PAHs that may have formed at early times by non-AGB sources. At later times, 
the effect of UV destruction will mainly effect the model by increasing the value 
of (wism)? which will now include the effective mass of the ISM gas that was 
cleared of PAHs by UV radiation, in addition to that cleared by shocks. 
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In summary, at early times the two processes discussed above cancel each 
other out, since the UV radiation can effectively destroy most of the PAHs 
that may have formed by interstellar processes before the first progenitors of 
AGB stars evolved off the main sequence. At later times, the two processes will 
be indistinguishable from those already considered in the paper: any interstellar 
PAH source will follow the time evolution of AGB stars, and any UV destruction 
mechansim will follow that of the destruction rate by SNRs. 

The global trend of PAH abundance with metallicity does not preclude the 
possibility of local variations of PAH abundance in individual galaxies. For 
example, the SMC which has on the average a metallicity that is 1/6 Z © (Dufour 
et al. 1982) has at least one region, the molecular cloud SMC Bl#l (Reach et al. 
2000), with a PAH-to-dust ratio that is comparable to the Galactic value (Li &; 
Draine 2002). Even in our Galaxy, there are regions with “super Galactic” PAH 
abundances (e.g. Ridderstad et al. 2006). Such local abundance variations are 
a natural consequences of processes, such as mixing of stellar ejecta and cycling 
between the ISM phases, that when globally averaged, will follow the general 
trend of PAH abundances with metallicity depicted in Fig. 18. 


4.4. The Evolutionary Trend of SN-Condensed Dust 

The far-IR dust content is in good agreement with the dust production by SN n, down 
to ~ 0.1 Z 0 (Fig. 18). Below this value, the model systematically overestimates the observed 
far-IR dust content. There may be several reasons for this discrepancy. 

First, we may have overestimated the gas mass, in the lowest metallicity sources. Indeed, 
the distribution of Hi of these galaxies extends farther out of the star forming region, as 
discussed in §4.1. We have attempted to correct this effect, for I Zw 18 and SBS 0335-052, 
by considering only the gas content associated with the optical galaxy. However, we could 
not correct for the gas located out of the star forming region and along the line of sight. 

Second, we may have underestimated the dust mass by overlooking a cold dust compo- 
nent. Indeed, Galliano et al. (2003, 2005) showed that the submillimetre excess observed in 
the SEDs of NGC 1140, NGC 1569, II Zw 40, and He 2-10, could be consistently explained 
by the presence of very cold dust, accounting for 40 to 80% of the total dust mass. On 
Fig. 18, the FIR-dust-to-gas mass ratio of these four galaxies is indeed, below the SN n 
production rate by a factor of ~ 2. Assuming that the high dumpiness of the ISM is a 
general property of low-metallicity systems, and that the filling factor and/or contrast den- 
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sity of these clumps rises when the metallicity drops, we have a natural explanation for this 
deviation. We can not address this issue, because of the lack of submillimetre data for these 
very low-metallicity galaxies. 

Third, the discrepancy between the predicted and observed SN II dust at low metallicity 
could be due to several parameters or assumptions used in the chemical evolution modeling. 
(1) The IMF could play an important role in the absolute value of the dust production by 
SN II. It would indicate that the slope of the IMF is metallicity-dependent, and that the 
contribution of massive stars is lower, at very low Z\su, which seems to be unlikely, both 
from a theoretical point of view (star formation), and an observational point of view (star 
counts). (2) A much more likely explanation could come from the fact that we assumed 
that the condensation into dust of the elements ejected by the SN II and their mixing in the 
ISM is instantaneous. If the major part of the dust was to condense into the ISM, and not 
directly into the SN Ii-ejecta, then the dust formation would be delayed after the death of 
the massive stars. (3) Finally, the dust production rates, computed from our dust evolution 
model, implicitely assume that the star formation rate of the galaxy is smooth and that the 
mixing is instantaneous. This hypothesis is certainly valid for evolved systems, but could be 
wrong for very young objects. For example, Legrand et al. (2000) suggested that the star 
formation history of I Zw 18 is not continuous. 


5. CONCLUSION AND SUMMARY 

The origin of the weakness of the mid-IR aromatic features - commonly attributed to 
emission from PAHs - in low-metallicity environments has been traditionally interpreted 
as the consequence of the increased selective destruction efficiency of the PAHs in these 
environments. In this paper, we presented a new interpretation for the observed correlation 
of the intensity of the mid-IR emission from PAHs with galaxies’ metallicity. In our model, 
the trend is a manifestation of the evolution of the abundance of PAHs that condensed in 
AGB stars. We also show that the inclusions of any source of PAHs, in addition 
to the AGB stars considered here, and a radiative destruction mechanism, in 
addition to the mechanical one by SNRs, will still retain the evolutionary trend 
of AGB condensed dust abundance with metallicity. 

To ascertain this trend and derive the abundance of PAHs, we modeled the UV to far-IR 
SEDs in a sample of 35 nearby galaxies, with metallicities ranging from 1/50 to ~ 3 Z Q . 
From these models, we determined the intrinsic stellar radiation field which was used to 
determine the abundance of the PAHs giving rise to the mid-IR aromatic features, and the 
abundance of the dust that gives rise to the far-IR. emission. We then used a chemical 



FIR dust 



Fig. 18. — Comparison between the metallicity trends of the PAH abundance derived from 
the observed SED and those derived from the chemical evolution model. The figure highlights 
the different evolutionary trend of SN ir- and AGB-condensed dust. 
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evolution model to calculate the abundances of SN- and AGB-condensed dust as a function 
of time - or metallicity. The model takes into account the delayed recycling of the ejecta 
from low-mass stars caused by their finite main-sequence lifetime. 

The main conclusions of this paper are the followings. 

1. In our analysis of the mid-IR spectra of nearby galaxies, we detected the 7.7 fim 
aromatic feature at the 4cr level in VII Zw 403, a dwarf galaxy with Z ~ 1/20 Zq. 
This is the lowest metallicity galaxy for which PAH emission has been detected to date. 

2. We find that the derived relative abundance of PAHs and FIR-dust evolves with metal- 
licity, suggesting distinct evolutionary histories for these two types of dust grains. 

3. Comparison of the observed dust-to-gas mass ratios with the results of our chemical 
evolution model shows the following: (1) the delayed injection of PAHs into the ISM 
provides a natural explanation for the paucity of these large molecules in low metallicity 
systems; (2) the subsequent rise in the PAH-to-gas mass ratio with metallicity is then 
a natural consequence of the increasing contribution of AGB stars to the chemical 
enrichment of the ISM as they evolve off the main sequence; (3) the evolutionary trend 
of FIR-dust-to-gas mass ratio with metallicity generally follows the distinct evolution 
of SN-condensed dust which is instantaneously injected into the ISM after the birth of 
the progenitor star. 

4. The inferred FIR-dust-to-gas mass ratio falls below the calculated value for SN-condense 
dust in the lowest metallicity galaxies. This discrepancy may be due to one or more of 
the following: (1) an overestimate of the Hi gas mass used to derive this dust-to-gas 
mass ratio; (2) an underestimate of the FIR-dust mass due to the possible presence of 
a cold dust component; (3) the gradual mixing of SN ejecta into the ambient medium; 
and finally (4) a more complex star formation history than used in the model calcula- 
tions. 

The success of our chemical evolution model, in reproducing the trend of PAH abun- 
dances with metallicity, strongly suggests the importance of stellar evolutionary effects in 
determining the abundances and composition of dust in galaxies. These will have important 
consequences for the opacity of galaxies and their reradiated thermal IR emission. Chemical 
evolution models for dust must therefore be an integral part of population synthesis models, 
providing a self-consistent link between the stellar and dust emission components of the SED 
of galaxies. 
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